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Abstract 

We use momentum twistors to evaluate planar loop integrals. Infrared divergences 
are regulated by the recently proposed AdS-inspired mass regulator. We show 
that two-loop amplitudes in M = 4 super Yang-Mills can be expanded in terms 
of basis integrals having twistor numerators. We argue that these integrals are 
considerably simpler compared to the ones conventionally used. Our case in point 
is the two-loop six-point MHV amplitude. We present analytical results for the 
remainder function in a kinematical limit, and find agreement with a recent Wilson 
loop computation. We also provide two-loop evidence that the logarithm of MHV 
amplitudes can be written in terms of simple twistor space integrals. 
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1 Introduction 



Recent years have seen tremendous progress in calculations of loop level amplitudes in 
planar M = 4 super Yang-Mills (SYM). State-of-the-art methods based on generalised uni- 
tarity [U [2] and sophisticated methods for the evaluation of loop integrals [3] have made 
the computation of various amplitudes at high numbers of loops and with many external 
legs possible, see e.g. [U El El Ej. Current techniques allow, for example, the numerical 
evaluation of six-point MHV or NMHV amplitudes at two loops [SJ [9]. Obtaining analyt- 
ical results for n = 6 or computing n > 6 amplitudes at two-loops however seems very 
involved using these methods. On the other hand, recent encouraging results [10] suggest 
that one should try to find analytical answers for amplitudes in M = 4 SYM, and despite 
the progress mentioned above it seems clear that for doing so we need new tools and insights. 

In [11] , an alternative regularisation for infrared divergences in planar M = 4 SYM was 
put forward (see [121 IS] f° r earlier references). It is inspired by the AdS description of 
scattering amplitudes [H] and consists in considering scattering amplitudes on the Coulomb 
branch of M = 4 SYM. The vacuum expectation values of some of the scalars give rise to 
masses that can be used to regulate the loop integrals. It was shown that this has a number 
of conceptual as well as practical advantages over dimensional regularisation, and that it 
facilitates the computation of amplitudes with many loops or external legs [151 US] • 

In this paper we propose a new method that combines the virtues of the mass regula- 
tor with those of momentum twistor variables |17^l . The latter are well-suited to describe 
planar four- dimensional scattering amplitudes. Recent papers discussed the one-loop box 
integrals using momentum twistors [T9~l [20] . Here we apply momentum twistors to higher 
loop integrals. We argue that integrals with certain momentum twistor numerators are 
much simpler compared to the integrals that are conventionally usecjE We show how to use 
twistor identities to express amplitudes quite generally in terms of our preferred basis. 

We focus on amplitudes with n > 6 external particles. The reason is that dual confor- 
mal symmetry, which is widely believed to apply to planar scattering amplitudes in M = 4 
SYM, entirely fixes the form of the four- and five-particle amplitudes, to any loop order 
[2~1~1 122] • (This is related to a conjectured relationship between Wilson loops and MHV 
scattering amplitudes [HI [231 [2U [21], for reviews see [251 [2S]). For n > 6 external legs 
there is the freedom of an arbitrary function, called remainder function, that can depend 
on conformal cross-ratios only. Indeed the two-loop MHV amplitude and corresponding 
Wilson loop have a non-trivial remainder function [23 El [2H]- As a specific example we 
then rewrite the six-point MHV amplitude in terms of our preferred basis integrals. We 
evaluate the latter and compute the remainder function analytically in special kinematical 
regimes. 

We are confident that the methods developed here and in our forthcoming paper [29] 
can be used to make contact with recent results for certain two-loop Wilson loops in various 

1 It is encouraging that the same variables were introduced independently at strong coupling [18] . 
2 We are particularly grateful to Nima Arkani-Hamed, Jake Bourjaily, Freddy Cachazo, Simon Caron- 
Huot and Jaroslav Trnka for sharing their ideas related to integrals constructed using momentum twistors. 
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kinematical regimes, see [30] and [3TJ [32] . For example, it is straightforward to apply the 
methods presented here to arbitrary n-point two-loop MHV amplitudes [33]. Moreover, 
while the correspondence with Wilson loops is, at the moment, restricted to MHV ampli- 
tudes, the methods we develop in this paper are quite general and can also be used to study 
non-MHV amplitudes. 

The outline of our paper is as follows: In section [21 we review the symmetries of scat- 
tering amplitudes focussing on momentum twistor variables and the massive regularisation. 
We give a one-loop example to show how twistor numerators can easily be dealt with. We 
then go on to describe in section [3] how to use twistor identities in order to transpose an 
amplitude expressed in terms of canonical integrals in terms of our preferred basis. We use 
the six-point two-loop MHV amplitude as an example, and discuss the virtues of the new 
basis. In section S] we evaluate the integrals for the six-point two- loop MHV amplitude and 
provide an analytical result for the remainder function in particular kinematical regimes. 
In section [5] we provide another example for the usefulness of twistor variables by showing 
that the logarithm of the two-loop four-point MHV amplitude can be written as a simple 
twistor integral. In section we present a simple but useful consistency test of loop inte- 
grands based in the soft limit. Finally, in the appendix we present more technical results 
quoted in the main text, as well as a proof that all leading singularities in M = 4 SYM are 
dual conformal covariant. 

2 Lightning review of momentum twistor space and 
massive loop integrals 

Dual conformal symmetry is an important property of planar colour-ordered amplitudes in 
M = 4 super Yang-Mills theory. Let us introduce the basic notions. Given the n incoming 
light-like momenta of a planar ordered amplitude, 

pf = \f\t (2.1) 
we define the dual coordinates in the usual manner [3U [33] , 

x aa_ x aa i=p aa (22) 

The dual Xi define a light-like polygon in the dual space. Under conformal transformations 
of the dual coordinates (dual conformal transformations) the polygon transforms but the 
edges remain light-like. This fact means that one can discuss dual conformal transfor- 
mations of planar ordered scattering amplitudes. In the M = 4 theory they are actually 
symmetries of the amplitudes, exact at tree-level [311 E3 EE] but anomalous when acting 
on a suitably defined finite part at loop level [2H [22J [35] . 

When taken together with the ordinary superconformal symmetry of the Af = 4 su- 
per Yang-Mills Lagrangian, the dual conformal transformations generate the Yangian of 
the superconformal algebra, Y(psl(4\4)) [39]. The Yangian is a symmetry of the tree- level 
amplitudes (with the usual caveats about contact terms in the action of the ordinary su- 
perconformal symmetry [101 SD S2])- At loop level the amplitudes break the ordinary 



4 



superconformal symmetry, however all leading singularities (see [I3j EH E5] for discussions 
of leading singularities) are ordinary superconformal invariants. This can be most easily 
seen in twistor space as they can be obtained by gluing tree-level amplitudes together in 
way which manifesty respects the superconformal symmetry [16J H?]. All leading singu- 
larities so far examined are also dual conformal. At one loop this essentially follows from 
the analysis in [IHJ |37] . We give a general proof of this fact for all leading singularities in 
appendix [A] This is relevant to the conjecture of [15] that all leading singularities can be 
obtained from an integral over a Grassmannian of a particular kind. This dovetails very 
nicely with the fact that, under mild assumptions, the Grassmannian integral can be shown 
to give the most general form of a Yangian invariant, i.e. that all Yangian invaraints are 
obtained from the integral with some choice of contour [^9l \50\ [51]. 

Now let us turn to dual conformal symmetry realised at the level of the loop integrals. 
In [35] it was noticed that the integrals appearing up to three-loops in the planar four-point 
amplitude formally exhibit a conformal symmetry in the dual space. Here we would like 
to discuss new integrals which also have this dual conformal property. Momentum twistors 
will be central to our discussion. These variables were introduced in [IT]. A very helpful 
discussion of the associated geometry is also given in [52J . They are the natural twistors as- 
sociated with the dual coordinate space which can be used to describe scattering amplitudes. 

A point in dual coordinate space corresponds to a (complex, projective) line in momen- 
tum twistor space. Two dual points are light-like separated if the corresponding lines in 
momentum twistor space intersect at some point in momentum twistor space. Thus the 
light-like polygon in dual space corresponds to a polygon in momentum twistor space with 
each line intersecting its two neighbouring lines as each dual point is light-like separated 
from its two neighbours. The n momentum twistors associated to this configuration of n 
light-like lines are defined via the incidence relations, 

%i = (^H i^i )) [A = x i \ct = X i+1 Xi a . (2-3) 

The momentum twistor transforms linearly under the action of dual conformal symmetry, 
as indicated by the fundamental sl(4) index A. Moreover the n momentum twistors de- 
scribing the polygon are free variables, in contrast to the dual points Xj which obey the 
constraints of light-like separation from their neighbours. The dual point Xi is associated 
the line described by the pair z\^_ i zf^ or (i — 1 i) AB for short. 

Note that only the lines (i — 1 i) correspond to the cusps of the dual polygon. One can 
consider other lines in momentum twistor space, for example the line (i i + 2). This line 
obviously intersects the four lines (i — li), + (i + li + 2) and (i + 2 i + 3) . It therefore 
corresponds to a point which is light-like separated from the four points x«, Xi + 2 and 
x i+3 . This distinguishes it as one of the two solutions to the four-particle cut conditions 
when a one-mass box integral is thought of in dual coordinate language. The other solution 
we will denote as + which is related to the first by parity. Both of the corresponding 
dual points are complex. 

The incidence relations allow one to express functions of the Xi in terms of momentum 
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Figure 1: Pentagon integral f|2.6[) with twistor numerator (13AB). The figure on the right shows 
the same integral in the more familiar dual space notation, with the dashed line denoting a 
numerator connecting to the complex point x a , see (|2.8p . 



twistors. For example we have 

2 (i — 1 i j — 1 j) 

x * = (i-uw-ij) ' (2 - 4) 

where the four-brackets and two-brackets are defined as follows, 

(ijkl) = e ABCD ZfzfZ^Z^ (ij) = XfXja . (2.5) 

The four-brackets are obviously dual conformal invariants while the two-brackets are in- 
variant under just the Lorentz and (dual) translation transformations. 

Since momentum twistors linearise the action of dual conformal symmetry it is clear 
that they are natural variables for discussing dual conformal integrals. Let us discuss 
the formulation of loop integrals using momentum twistors. Recent papers have already 
employed these variables to discuss the one- loop box integrals [191 12H] • Here we will describe 
how other integrals can be similarly formulated using momentum twistors. In particular 
we are interested in integrals which make explicit use of the twistor lines which do not 
correspond to the dual X{ like the line (i i + 2) described above. As an example we will 
consider a one-loop pentagon integral, 

d 4 Z AB {ABU) (4512) (2345) 

(2.6) 



(AB12)(AB23)(AB34)(AB4:5)(AB51) ' 



Here the integration is over the space of lines in momentum twistor space (AB). As we 
have discussed, this is equivalent to an integration over points in dual space. Recall that 
the momentum twistors are only denned up to scaling so it is important that the above 
expression has zero scaling weight in all of the momentum twistor variables. Taking into 
account that the integration measure 

d A Z AB ~d A x (AB)\ (2.7) 

has scaling weight 4 for the twistors A and B, we see that (12.61) indeed has zero scaling 
weight. The numerator of the integrand in ( 12. 6ft contains the factor (AB13) depending on 
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the line (13) which, as we have discussed corresponds to a complex point. 



The integral form in (12. 6p is manifestly dual conformal invariant, being written in terms 
of four-brackets, even though we are describing a pentagon integral. In dual coordinate 
notation the pentagon integral is 



x 52 x 35 f d A Xp /„ 

9 2 2 2 2 2' V 



x 5a J Ln X 01 J '02 J '03 J '04 J '05 

where the point x a is one of the two solutions to 

X 5a = X la = X 2a = X 3a = ■ (2-9) 

We can see from (12. 8p that the integral form is dual conformal invariant, having dual con- 
formal weight zero at all of the points X\, x%, x^, X4, x$, x a . The existence of a sixth point x a 
was crucial in being able to make a dual conformal pentagon integral. It is not independent 
of the other five dual points Xi, . . . ,x$, being defined as a solution to the (one-mass box) 
cut-conditions (12.91) . nor is it real. The integral (I2.6P is therefore not invariant under parity, 
i.e. it is chiral. 

The integral (12.61) naturally leads us to introduce another tool that we will need in this 
paper. It is infrared divergent so we will need a regulator. We would like to preserve the 
dual conformal symmetry so we will use the AdS regular isat ion introduced in [UJ. For 
actual calculations we will use the regularisation where all masses are equal. In particular 
this means the outermost propagators in the planar loop integrals that we are studying are 
modified as follows, 

4-— > 2 I 2 - (2-10) 

X ij X ij + 171 

In our example of the dual conformal one-loop pentagon this means that the five prop- 
agators need to be altered to incorporate the mass terms. In momentum twistor language 
this can be achieved by turning each of the pairs of twistors associated with the dual Xj 
into a bitwistor with a component proportional to the infinity twistor I AB . Thus we make 
the replacements in the denominator of (12. 6p 

Zf_ x Zf -> Xt% = Zf_ x Zf + m\ % - 1 , (2.11) 

where [•, •] stands for the antisymmetric part. This has the effect that each of the propagator 
factors becomes 

(2.12) 



(ABi-li) (ABi-li) (ABX^u) (ABi - 1 i) + m 2 (AB)(i - 1 i) 

We will leave the numerator of (12. 6 j) unchanged. In fact the definition of the numerator 
factor does not matter in this case as any potential modification would be 0{m 2 ) and the 
integral is only logarithmically divergent so these terms can be safely ignored for the present 
purpose. Thus we arrive at the final definition for the regularised pentagon integral, 

d 4 Z AB {AB\$) (4512) (2345) 

i% 2 (AB12)(AB23)(AB3A)(ABA5)(AB51)' [ ' } 
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Let us proceed to find a Feynman parametrisation of this integral. In order to do this it is 
helpful to normalise the denominator factors so that they carry no scaling weights in the 
external twistor variables. We will define 

yAB y\A 7 B] 

yAB _ _ , 2 tAB f -i *\ 

Then the integral can be written 

(4512) (2345) f d 4 Z AB (AB13) 



Ili<ii + 1) 7 ^ 2 (A j BXi3)(ABX 23 )(ABX34)(ABX4b)(ABX 51 ) 
Introducing Feynman parameters we can write this as 

_ 4! (4512) (2345) f d*Z AB f n4 (AB13) 



(2.15) 



+ J in 2 J D a (ABY) 5 ' (2 - 16) 

where the bitwistor Y is defined by 

Y = a 12 X 12 + a 23 X 23 + a u X u + a 45 X i5 + a 51 X 51 . (2.17) 

We write the measure of the Feynman parameter integration as D 4 a to remind ourselves 
that this is a projective integral. The numerator factor can now neatly be written as a 
derivative, 

3!(4512)(2345) f f *Z„ 1 



+ 7 ' J in 2 (ABY) 4 ' 

The integral over the line (AS) can now be performed as described in [191 120] with the 
result 

, (4512) (2345) f 1 (4512) (2345) / (13F) 

+ 7 ( YY r U i (^ + 1 ) J ( YY r 

The remaining four-brackets in the numerator and denominator can now be expressed in 
terms of the twistor variables, 



(13Y) = a 45 



(1345) , 2 



+ m 2 (13) 



(2.20) 



(45) 

(YY) = 2 ®i~ii®i+n+2xli +2 + ^2 a »+i m2 • ( 2 - 21 ) 

i i 

Finally combining everything and ignoring the 0(m 2 ) term from the numerator in f)2.20p 
we obtain 

2 2 2 f d^a5{J2i^n+i - l)«45 „, 2^ 



J 5 = -2^*^35 / — ^ — '- " + 0{m*) . (2.22) 

J ai-ua i+1 i+ 2 xl i+2 + m 2 \ 

So in summary we see that the twistor numerators can be easily dealt with. In the next sec- 
tion we are going to explain how such integrals can be used in order to simplify expressions 
for loop integrands. 
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^6;2;1 h;2;5 h;2;9 h;2;12 

Figure 2: Representative integrals appearing in two-loop MHV amplitudes. In the pictures the 
specific choice n = 6 was made. The dashed lines stand for numerator factors that depend on the 
loop momentum. 



3 New representations of loop integrands 

In the previous sections we reviewed how momentum twistors and a massive infrared reg- 
ulator can be used to discuss planar loop integrals. We also showed that certain non-local 
twistor numerators can easily be treated in an approach using Feynman parameters. We 
will now explain how to use these non-standard numerators in order to derive simpler ex- 
pressions for the integrand of multi-loop and multi-leg amplitudes. 

As a case in point we are going to discuss the two-loop six-point MHV amplitude. The 
form of the two-loop six-point MHV amplitude given by Bern et al [8] and transposed to 
the massive regularisation is 



6 16 



1 13 

77 J2 ^Q/ 6 ;2; 4 + 0(m 2 ), (3.23) 



12 perms/flips 



with the integrals correspond to the /W of [8j, except that we include a normalisa- 
tion factor in order to make them dimensionless, and the coefficients Cj, given in [Sj, are 
modified accordingly. |f| The integrals appearing in ( I3.23P are of the double box, pentabox 
and double pentagon type, representatives of which are shown in Figure |5J The latter two 
involve certain numerator factors that depend on the loop momentum. Since they are more 
complicated than the double box integrals, we wish to eliminate them in favour of simpler 
integrals. We will now explain how this can be done quite generally. 

Consider for example the following six-point pentabox integral, 

f dZ AB dZ CD (AB12) 
6;2;9 K J (AB23) {ABU) (AB45) (AB56) (ABCD) (CD56) (CD12) (CD23) ' ^' 2 ^ 

which is shown in Fig. |2j We would like to exchange the numerator (AB 12} for terms 
which either cancel a propagator or give a preferred 'non-local' type numerator of the form 
(AB35) or (AB35). We can do this by writing 

(12) = 6 23 (23) + 634(34) + 6 45 (45) + M56) + 6 35 (35) + 635(35) . (3.25) 



3 In [8] there are also two further integrals depending explicitly on /1, the (— 2e)-dimensional component(s) 
of the loop momenta in dimensional regularisation. It has been observed that these terms cancel to 0(e) in 
the log Mq to the two-loop level, which suggests that analogous terms, if present in the mass regularisation, 
should be 0(ra 2 ) [16]. 
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Figure 3: Example of identity (|3.25p when used to expand /6;2;9 in terms of our preferred basis. We 
do note display the (loop-momentum independent) prefactors and normalisations of the integrals. 



Here 

(35) AB = (234 ■ ) A ( ■ 456) B - (A, B) . (3.26) 
By projecting with different twistors we can solve for the coefficients. We obtain 

(1245) (6245) (5312) (6235) (1234) 

23 ~ (2345) ' 34 ~ (6345) (2345)' 45 ~ ~ (6345) (2345) ' 

(1234) (6245) (1234) (1235) 

56 (6345)' 35 (4635) (2345)' 35 (6345) (2345) ' 1 ' 

Thus we have expressed the pentabox integral / in terms of double box integrals and 
pentabox integrals with the preferred numerators (AB35) and (AB35), up to +0(m 2 ) 
terms that arise due to the difference between (. . .) and (...). Moreover it turns out that 
the integral with (AB35) in the numerator is equivalent, again to 0(m 2 ), to the (AB35) 
integral. The resulting identity is schematically shown in Figure [3j 

One can use ( 13.251) and (13.271) and their reflected versions and similar identities obtained 
upon rotation to expand any unpleasant numerator. For example if we consider the double 
pentagon integral with numerator (AB12) (CD3A) (with AB and CD representing the loop 
integration variables) we can expand it by using ( I3.25P for the first factor and the similar 
identity obtained by reflection and rotation for the second factor. The above are examples 
of a completely general identity which can be used for any such integral. It can be found 
in the Appendix, and can be used to simplify any of the pentaboxes or double pentagons 
appearing in any two-loop amplitude. 

There is another type of term in the two-loop amplitudes that we would like to remove, 
namely the "kissing box" topology, which are a product of two one-loop box integrals, see 
l6;2;i in figure [2j Of course these are analytically quite simple, but we observe that they 
come with certain prefactors such that they contain rational factors, a feature which we 
would like to eliminate. We can rewrite those integrals as a double-pentagon integral by 
multiplying their integrand with 1 = (ABC D) / (ABC D) , and use the following identity to 
decompose the (ABCD) in the numerator: 

(ikAB) (jlCD) + (jlAB) (ikCD) = (i jAB) (klCD) + (klAB) (i jCD) 

-(jkAB)(UCD) - (UAB)(jkCD) - (ijkl)(ABCD) . (3.28) 
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h;2;ml h;2;m2 h;2;m3 -^6;2;7 ^6;2;15 

Figure 4: Integrals appearing in the improved representation of the six-point two-loop MHV 
amplitude, c.f. (13.291) . The dashed lines indicate non-local twistor numerators, and a normalisation 
factor is not shown in the figure. Both are given explicitly in the Appendix. 



The identity holds for any i,j,k,l. The integrals we obtain in this way can in turn be 
reduced using (I3.25p . as described above. 

The numerator identities described above can be used to obtain an improved represen- 
tation of the two-loop six-point MHV amplitude, which is given by [53J 



4 

12 perms/flips 



^6;2;ml — —h;2;m2 ~ T-^6;2;ra3 + T-^6;2;7 + 2^6;2;15 



+ 0(m 2 ), (3.29) 



where all integrals are depicted in Figure HI and their explicit definitions including their 
normalisation is given in the Appendix. The new representation ( 13 .29 j) is manifestly more 
compact than (I3.23p . as it requires only 5 as opposed to 13 integrals. Moreover, it has a 
number of other virtues, as we discuss presently: 

• We find that all integrals appearing in ( 13 .291) . to the order that we evaluated them 
analytically, do not contain any rational functions multiplying the (poly-) logarithms 
appearing in them. Experience shows that this is a generic property of scattering 
amplitudes in M = 4 SYM, and it is remarkable that it holds for all integrals in our 
basis, even before summing them up. This is not the case for the integrals appearing 
in f)3.23p . see e.g. the expression given for i6;2;ii i>e. I^ 1 ' in Appendix A. 2 of [8]. 

• The previous point hints at a hidden simplicity of the integrals appearing in ( 13 .29j) . 
This can be explained in part by a new method based on differential equations 
that will be presented in a forthcoming paper 



• It has become standard to evaluate loop integrals using Mellin-Barnes representations 
[3J. We remark that the bottleneck of the calculation of n > 6 two-loop amplitudes 
used to be the double pentagon integrals, as e.g. /6;2 ; i2 and /6;2 ; i3 in f l3.23[) . This 
becomes clear from inspecting the dimensionality of the Mellin-Barnes representations 
used to evaluate these integrals: the ones used for those two integrals in [8J are 18- 
fold. In the improved representation (I3.29p . we need 12- fold Mellin-Barnes integrals 
at most0 

Two remarks are in order: Firstly, the method we described here to simplify the repre- 
sentations for loop amplitudes/integrands applies more generally. For example, it can be 



4 The dimensionality of the Mellin-Barnes integrals is further reduced considerably when taking the 
regulator limit, i.e. e -> in dimensional regularisation, or m 2 — > here. 
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used to simplify the expressions for n-point two-loop MHV integrand [33] , the recently pre- 
sented form of the six-point two-loop NMHV integrand [9], or any other loop amplitudes. 
Secondly, we have already argued that the expressions we gave for the six-point two-loop 
MHV integrand greatly facilitates its evaluation, and we will see in the next section that 
this is indeed the case. However it may very well be that ultimately there exist yet better 
representations (from the point of view of evaluating the integrals), that can be obtained 
for example by replacing the remaining double box integrals Iq-2-,7 and /6;2;is by pentabox 
integrals with "magic" numerators (other than the ones already appearing in I^mi and 
h;2-,m3)- We leave this question for future work. 

4 Two-loop amplitudes in J\f = 4 SYM 

Here we discuss the evaluation of the loop integrals in the new basis described in the 
previous section. As a specific example, we will present analytical results for the six-point 
two-loop amplitude in certain kinematical limits. We stress that the methods presented 
here can be straightforwardly used to evaluate two-loop amplitudes with more external legs 
or different helicity configurations. 



4.1 Structure of the loop corrections for MHV amplitudes 

Let us begin by recalling the expression for the one-loop six-point amplitude. Its integral 
representation is pQ 



M, 



(i) 



1 E 

6 perms 



2 



+ 0{m 2 



(4.30) 



and using the explicit expressions for the integrals in equations (1C.94j) and (1C.95j) we see 
that it is given by (dropping 0(m 2 ) terms) 



M. 



(i) 



1 6 



i=l 



log 2 



m 



x 



i,i+2 



+ F, 



(i) 



7T" 

T + 2 



1 6 



i=l 



iy 2 ™ 2 

- log ~i — log 



1 



X i,i+3 



X 



+ 7 l°g 2 

4 xj 



2 X i,i+3 



i,i+3 "°i+l,t+4 
2 



Li 2 

2 



(4.31) 



1 - 



, y» 2 ,y. 2 

• x 'j,j+2- x 'j+3,j+5 
/v»2 ,-yi 2 

x V+3 x i+2,i+5 / J 



where we used the identity Li 2 (l — x) + Li 2 (l — 1/x) + 1/2 log x = 0. We remind the reader 
that the function satisfies the anomalous Ward identities of 



We now wish to investigate this amplitude at the two-loop order. It is worthwhile to 
discuss the expected structure of n-point MHV amplitudes at higher loops. In reference 
[16J it was suggested that the logarithm of the loop corrections to MHV amplitudes in the 
massive regularisation should have the following form, 



logM n 



E 



8=1 



7(fl) log 2 



16 



x 



i,i+2 

m 2 



Go(a) 



log 



x. 



iji+2 

m 2 



+ f(a) 



+- 7 (a)F« +1Z n ( Ui ,a) 



C{a) + 0(m 2 



(4.32) 
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This is motivated by the ABDK/BDS ansatz [HE], allowing for a function lZ n of conformal 
cross ratios. Let us briefly discuss this formula. Firstly, it states that the infrared divergent 
part of the amplitudes, which manifests itself through logarithms in m 2 as m 2 — > 0, should 
exponentiate, and that this exponentiation is governed by certain anomalous dimensions^] 
7(a) and Qo(a), in accord with [56j |57]. The latter take the values 



7 (a) = 4a - 4( 2 a 2 + 0(a 3 



Go(a) = -Csa 2 + 0(a 3 



(4.33) 



Secondly, the finite terms in (I4.32p can be thought of as an explicit solution to the Ward 
identities of [22], with ^7(0) Fn being a particular solution. The full finite part allowed by 
the Ward identities is the particular solution plus a function of the available conformal cross- 
ratios TZ n (ui). It vanishes for n = 4 and n = 5, since non-vanishing conformal cross-ratios 
appear only starting at n = 6. Thirdly, the specific choice of the kinematic-independent 
terms f(a) and C(a), determined by the four- and five-point case, was made in such a way 
that the remainder function vanishes in the collinear limit. They are given by [16] 



c 



f(a) = ^a 2 + 0{a A ) , 



C( a ) = -^a 2 + 0(a 3 ) 



(4.34) 



Specialising to the six-point case at two loops, we see that Fq 1 ^ was defined in the second 
line of (I4.3ip and TZe(a) = 0(a 2 ) by definition. As we discussed, equation (I4.32p implies 
that all logarithms in m 2 of are captured by 



-M, 



(i) 



1 1 2 

log 

16 & 



x 



i,i+2 

m 2 



(2) 



log 



•^1,1+2 

m 2 



(4.35) 



Moreover, we have that 



nt\u x ,u^) = M 6 (2) - Eqn. gSSJ) - V> F 6 (1) - ^-(4 , 



12 



(4.36) 



where 



111 



™,2 ,-y»2 

x 13 x 46 

2 2 ' 
x 14 x 36 



U 2 



,7^2 ,-y»2 

x 24 x 15 

2 2 ' 
x 25 x 14 



/y>2 ,-y»2 

x 25 x 26 

2 2 ' 
x 36 x 25 



(4.37) 



are the three non-vanishing cross-ratios at six points. In the following section, we will 
compute the loop integrals contributing to M\ and by virtue of (I4.36P obtain a result for 



K 



(2) 



4.2 Evaluation of two-loop integrals and amplitudes 

Let us explain the general strategy for computing the integrals contributing to the two- 
loop amplitude. As was explained in section |3l the twistor numerators can be conveniently 
dealt with, and the introduction of Feynman parameters for the propagators in the usual 
way allows to carry out the loop integration. It is then straightforward to write down 

5 The notation 7(a) = 2r cusp (a) is also widely used in the literature to denote the cusp anomalous 
dimension 1551 . 
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Mellin-Barnes representations for all integrals. The m 2 — > limit can be implemented by 
MBasymptotics [58], and barnesroutines [5H] can be used to further simplify the result- 
ing expressions by using Barnes lemmas where applicable. The resulting expressions are of 
similar complexity as those obtained in the Wilson loop computation of [60] (recall that 
MHV amplitudes are conjectured to be dual to n-cusped light-like Wilson loops). One 
may therefore hope to derive analytical answers for the remainder function, for example 
by exploiting that the cross-ratios that the latter depends on are invariant under certain 
Regge limits, as was successfully done in [6U] • A full analytical evaluation of the remainder 
would go beyond the scope of this paper. Here, we content ourselves to provide analytical 
results for the six-point remainder function in certain kinematical limits^. 

The Mellin-Barnes representations we obtain can also be readily evaluated numerically. 
It is amusing to note that the double-pentagon integrals of [SJ , which constitute the bottle- 
neck of the calculations in the conventional basis, have been replaced by simpler penta-box 
and double pentagon integrals, and the finite double pentagon integral with "magic" nu- 
merator, which is much better behaved (In fact the latter integral is infrared finite. It would 
not be an exaggeration to say that the momentum twistors allow us to replace the most 
difficult integral in the conventional representation by the simplest integral in the new rep- 
resentation). Therefore we are in a position to produce high-accuracy numerical results for 
the remainder function at arbitrary (Euclidean) values of the kinematical variables, using 
e.g. the MB package [6Tj . We remark that it is straightforward to extend this statement to 
other two-loop amplitudes, e.g. of higher multiplicity or for other helicity configurations, 
since the contributing loop integrals are expected to be of the same type. Therefore in 
particular we have also numerical access to non-MHV amplitudes^, which was not possible 
previously as the duality with Wilson loops is restricted to MHV amplitudes so far. 

Let us now return to our specific six-point MHV example. In section [3] we have already 
alluded to the fact that apart from containing fewer contributing integrals, the new repre- 
sentation of the loop amplitudes also has other virtues. In order to see this, let us evaluate 
the integrals in the improved representation of section |3j To begin with, we verified the 
infrared structure (I4.35p . The reader is invited to check this by the explicit formula we 
give for the log 4 m 2 and log 3 m 2 terms of the integrals in the appendix. We also computed 
analytically the log 2 m 2 terms but refrain from writing out the lengthier formulas. 

For simplicity, let us now consider the symmetric kinematical configuration where xj i+2 = 
1 and xj i+?j = w -1 / 2 , which sets all cross-ratios equal Ui = u. As an illustration of the above 
statements about the infrared structure of the amplitude, let us quote the result for the 

6 A new method for obtaining analytical results for loop integrals based on differential equations [33 EI] 
will be presented in a forthcoming paper [29]. 

7 Preliminary results for the six-point NMHV amplitude in dimensional regularisation where presented 
in 0. 
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kinematical point u = 1, for which we obtain (denoting L = logm 2 j^| 



h;2;ml = -6& + 0(m 2 ), (4.38) 

h;2;m2 = ~ ^ ~ ^ ~ ^ L ~ 1 ^4 + 0(m 2 ) , (4.39) 

h;2;m3 = 0(m 2 ) , (4.40) 

J 6 . 2;16 = ^L 4 - 4( 2 L 2 - 16( 3 L - 12C 4 + 0(m 2 ) , (4.41) 

J 6;2;7 = ^L 4 + 2( 2 L 2 + 12(3-^ + 12C 4 + C(m 2 ) . (4.42) 
At the symmetric point ( I3.29P becomes 

M 6 (2) = 1 [3/ 6;2;ml - 6J 6;2;m2 - 3/ 6;2;m3 + 3/ 6;2;7 + 6/ 6;2;15 ] , (4.43) 

Q 7 97 

= -L 4 — -L 2 — 3(3-^ + — — 7r 4 + 0(m 2 ) . (4.44) 

8 8 S 160 V y K 1 

and we find perfect agreement with the expected infrared structure according to equation 
(I4.35p . Moreover, from (I4.36P we see that 

W? ) (l,l,l) = — , (4-45) 



in perfect agreement with the corresponding analytical result obtained in [60] for the dual 
Wilson loop. We also evaluate the limits where u — > and u — > 00 analytically, and obtain 
(the intermediate results are given in the Appendix) 

7T 2 177T 4 

hm Tlf\u, u,u) = — log 2 u + — - + 0(u) , (4.46) 
u->o 1440 



and 

4 

hm ^ 2) (m,m, M ) = ~ + , (4.47) 

moo 144 

again in perfect agreement with the Wilson loop result of [60] . 

In summary, the new representation of the loop integrand, combined with the virtues of 
the regulator of [TT] . has allowed us to easily produce an analytical result for the six-point 
remainder function (in a kinematical limit). It seems straightforward to evaluate other 
two-loop amplitudes in a similar way. 

5 Twistor space integral for exponentiated four-point 
two-loop amplitude 

In the previous sections we argued that momentum twistors can be used to simplify the 
integral representations for two-loop amplitudes. In particular, we exploited a new form 



8 All coefficients were obtained analytically, except the constant in Ie-,2;mi, which was obtained to high 
numerical accuracy. 
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of the two-loop six-point MHV amplitude to evaluate the remainder function and obtained 
analytical results in certain kinematical limits. 



Here we wish to focus on MHV amplitudes and give a further application of momentum 
twistors. From the general structure of infrared divergences, see equation ( I4.32[) , it is clear 
that a good quantity to study is the logarithm of the loop corrections M n of MHV ampli- 
tudes. We recall that the mass regulator as such is already convenient for this purpose, as it 
allows to switch between M n and log M n without having to keep 0(m 2 ) terms^l However, 
the L-loop amplitude starts with log 2i m 2 as most divergent term, whereas in log M n 
the most divergent term is log m 2 , and so many terms must cancel between different loop 
orders when taking the logarithm. One can ask whether an integral representation exists 
that directly yields this logarithm, without having to compute those higher terms in the 
first place? 

As we show presently for the four-particle amplitude at two loops, this is indeed the 
case. There exists a natural way to present logM4 that has the above properties. As we 
will see, momentum twistors will again be crucial. 



To begin with, we recall the expression for the four-point amplitude to two loops, 

(5.4? 



(logM 4 ) (2) = i 



I 2 (s,t,m 2 ) + I 2 (t,s,m 2 ) - - (/i(s,t,m 2 )) 2 



where Ji and l 2 are the one- and two-loop ladder integrals, respectively, and s = xf 3 and 
t = X24. The result for the small m 2 expansion is in agreement with (I4.32p for n — 4, and 
is explicitly given by [11] 

(log M*)« = C2 log ^ log ^ - Cs (log ^ + log ^ - Att 4 + 0(m 2 ) . (5.49) 

The idea is now to combine all integrals on the r.h.s. of ( 1548ft and to simplify the resulting 
expression using momentum twistors. At this point we recall an identity that was given in 
section El 

(1345) (2ACD) + (2445) (13O0) = (1245) (34C5) + (3445) (12CD) 

-(2345)(41C5) - (4145)(23C5) - (1234)(450D) . (5.50) 

We see that, after symmetrising the integrand in the integration variables the coefficients 
of the terms exactly such that we can use (I5.50p ! In doing so, we have to take into account 
that all of the propagators in the integrals except the middle line of the double ladders 
contain a +m 2 term. Let us introduce the notation (1234) = (1234) + m 2 (12)(34) and 



dZ AB dZ CD Q (1234) 



3 



I[Q] ■= I 77-™ " * - : ,. ' ^7TT^ ■ (5-51) 



[ITT 



YT i= i [<i — 1, iAB)(i - 1, iCD)) (ABCD) 



where the cyclicity of the indices is understood and Q can be a loop-dependent numerator 
that is to be included in the integration. Then, we can write 

(logM 4 ) (2) = -I[(13AB)(24CD) + m 2 (1234)(45)(C5) + 45 ^ CD] . (5.52) 



3 In dimensional regularisation, 0(e) terms are required for this. 
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This is the main result of this section. The logarithm of the four-particle amplitude at 
two loops can essentially be written as a single twistor integral, with "magic" numerator 
(13AB) (2ACD). In what follows we will verify explicitly that (I5.52p gives the correct an- 
swer (I5.49p . As advertised, we will not have to compute any terms higher than log 2 m 2 . 

A word of caution is due with respect to the m 2 (1234) (AB) (CD) term in the numerator 
- usually the integrals we deal with have only logarithmic divergences in m? as m 2 — > 0, 
so that such terms remain 0(m?) after integration and can be neglected. Here however, 
writing all terms over a common denominator has increased the number of propagators and 
lead to a more singular denominator. As a result this term has to be kept, as we will see. 

We now compute the small m 2 expansion of the two integrals in (15.521) . Proceeding as ex- 
plained in section[2]we find the following Feynman parameter formula for I 2 [(13AB) (2ACD)], 

2+2 r°° cW(ELiA-i) r^(E-=i«i-i)«5 , 2 , 



23 f Jo Wss + + m 2 ) J ( 7l73S + W + ™ 2 ) 3 + ° {m } ' (5 ' 53) 

where ji = ci^fti + The small m 2 limit is then readily evaluated using Mellin-Barnes 
methods, with the (remarkably simple) result 

2 2 / 2 2 

m , m , , „ ^ , / , m , m 



4C 2 log — log — + (-4C 3 + 8( 2 ) log — + log — - 17.768 + 0(m 2 ) . (5.54) 
s t \ s t J 

The m 2 (1234) (AB) (C D) term can be evaluated in a similar way. We find 

log— + log — j - 11.454 + C(m 2 ). (5.55) 

Combining equations f !5.54p and f)5.55p according to (I5.52p we have complete agreement 
with equation (I5.49p . It is remarkable that no log 4 m 2 terms appeared in the calculation, 
and no cancellations of say Li4 terms were required. 

In conclusion, we found a remarkably simple integral representation of logM 4 at the 
two-loop level. The fact that we were able to use the twistor identity ( I3.28P relied on the 
relative coefficient between the one-loop and two-loop integrals, which in turn is fixed by 
the infrared structure of the amplitude. This may suggest that similar formulas can be 
written at higher loops or for higher multiplicities. 



6 Consistency of the loop integrand with soft limits 

It is well known that collinear and soft limits can provide an important consistency check 
for loop amplitudes pp. Here we propose that the soft limit can also be used in a more 
direct way as a consistency check of the loop integrand. This check works in a very simple, 
almost pictorial way, and does not require any loop integrals to be evaluated. 

How do the scattering amplitudes behave in the soft limit, where where the momentum 
of one of the scattered particles is sent to zero, e.g. p n — > 0? The generic structure is, e.g. 
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at the one- loop level pQ, 

lim M n M n _! + div , (6.56) 

where the div is some universal soft divergence We argue that the divergent term in (I6.56P 
arises due to the non-commutativity of two limits: the soft limit and the regulator limit 
(e.g. e — > in dimensional or m? — > in mass regularisation). For finite value of the 
regulator, there is no reason for a divergent term to appear, and we expect the soft limit 
to be completely smooth, 

lim M n -> M n _x , (6.57) 

Pn->0 

to all loop orders. 

How will we be able to use (I6.57P in practice? We would like to argue that (I6.57P can be 
a useful guiding principle for constructing the correct loop integrand. Indeed if we can find 
an ansatz for the amplitude that manifestly respects 06.571) at the level of the integrand, 
then it will trivially also be true for the integrated amplitude. 

Let us now give some examples to see the usefulness of this approach. Our first example 
is the one- loop MHV amplitude M n . Let us consider the ansatz 

M n = ~ ^ F 2me(Xa, X a+ l,X b , X b+1 ) , (6.58) 
a,b 

where the sum runs over all inequivalent unordered pairs (a,b). This is indeed the correct 
answer [lj, at least up to +0(m 2 ) terms, which may in particular contain parity odd 
contributions. Note that some of the "two-mass easy" integrals in ( I6.58P are in fact "one- 
mass" integrals (for specific values of the labels), but we shall leave this distinction implicit. 
Let us now analyse the consistency of (I6.58P with the soft limit ( I6.57P . The integrand of a 
generic two-mass easy function F 2me is given by 



! ,.2 ,.2 ,.2 ,.2 

x a+l 



x a+1,6+1 J/ b+l,a x a+l,& 



X h+ i < y a0 ~o+i,a~a+i,o ^ (6.59) 

+ m2 ){xl+l,i + m2 )( X b,i + m<2 ) ( X b+l,i + m 2 ) ' ' 



X n 



where Xi is the integration variable. We see that there are two different possibilities: in 
the limit where a momentum belonging to a massive corner vanishes, we simply obtain 
lower-point integral i^me, with the unchanged labels. On the other hand, if the momentum 
taken to zero is p a or pb, the numerator in (16.591) makes the expression vanish. This is 
important, because otherwise we would have obtained an unwanted triangle integral with 
doubled propagator. So we see that one-loop MHV amplitudes (I6.58P manifestly satisfy 
d637D. 

Our second example is the two-loop MHV amplitude, which works in exactly the same 
way. Since the improved representation uses twistor variables, we have to formulate the 
soft limit in terms of those variables. When momentum p n goes to zero, we have that the 
twistor Z n becomes a linear combination of Zi and Z n _i, i.e. 

Z n aZt + f3Z n _! . (6.60) 
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The interested reader is invited to verify that the integrand of the five-particle amplitude 
can be recovered in this way. 

We remark that the way the MHV amplitudes manage to be manifestly consistent with 
( I6.57P is very reminiscent of how the same soft limit works for tree-level amplitudes, as was 
discussed in [38]. So far we only discussed the soft limit that does not change the helicity 
class of the amplitude. It is expected that the soft limit that goes from an n-point N fc MHV 
amplitude to an (n — l)-point N fc_1 MHV amplitude is also smooth. This should provide a 
useful consistency check on the loop integrand of non-MHV amplitudes. 

We hope to have convinced the reader of the usefulness of this test. A word of caution 
is due: it is well-known that the soft/collinear limit cannot completely constrain the ampli- 
tude, as there can be non-trivial functions that vanish in all limits, the remainder function 
TZq being a good example. Similarly, one can easily think of integrands that vanish in all 
soft limits due to numerator factors, such as e.g. (1234) (3456) (5612) in the six-point case. 
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A Dual conformal symmetry of leading singularities 

Here we give a proof that all leading singularities in planar M = 4 super Yang-Mills theory 
are dual conformal covariant. The calculation proceeds as follows. First we will start with 
an arbitrary planar graph constructed from tree-level amplitudes glued over the internal 
edges. If there are 41 cut propagators then the fact that the internal lines are on shell cut 
propagators means that the loop integration is localised by these conditions to a discrete 
set of possibilities of which one must be chosen. If it there are fewer than 41 cut propagators 
then one must specify further the nature of the contour - there will be other poles around 
which the loop integration must be taken - if not then the leading singularity would be zero. 
Such leading singularities are called composite leading singularities. Once the contour has 
been completely specified as a closed contour integral all internal momenta are fixed in 
terms of the external ones. 

The vertices of the graph correspond to tree-level superamplitudes of M = 4 super Yang- 
Mills theory. These tree-level amplitudes can be of any helicity type (MHV, NMHV, etc) 
including the three-point MHV amplitude. Following [46] we will replace any non-MHV 
amplitude with its BCF expression, e.g. for an NMHV tree-level amplitude we will write 
it as a sum over the gluings of four tree-level vertices of lower MHV degree. In this way we 
will increase the number of loops but will reduce every graph to a sum over graphs with 
only MHV and MHV3 vertices. These graphs are called primitive in the language of |46j . 
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Now that the graphs have been reduced to ones involving only MHV and MHV3 vertices, 
the next step is to perform the sum over exchanged states on each internal leg. In other 
words we have to perform the Grassmann integral associated to each internal leg. The 
Grassmann integrations can be performed via three different operations which we will call 
A,B and C: 

A. MHV 3 glued to MHV: 

The relevant Grassmann integral is 

J d%5\ Vl [12] + Vl [2l] + V2[H])5 8 (-XiVi + q K ) = [12] 4 5 8 (g). (A.61) 

In other words as far as the Grassmann factors are concerned one can replace the MHV3- 
MHV n _i pair with a single MHV n vertex. 

B. MHV glued to MHV: 

Here the relevant Grassmann integral is 

J d\5 8 (q Kl + X lV i)S 8 (-XiVl + Qk 2 ) = 5 8 (q)5 4 ((lq Kl )) = 5 8 (q)5\(lq K2 )) (A.62) 

Thus one obtains an MHV vertex with a decoration factor of the form ^((/q^)). Under 
the delta function S 8 (q) we can introduce a dual notation by writing 

5*i =ipi-^2 = -Qk 2 - (A. 63) 

One of the tp variables will be identified with a external dual 6 variable if one or other 
side of the gluing line corresponds to the outside of the original graph. Otherwise they are 
just shorthand for a given linear combination of supercharges, e.g. qx 1 which may contain 
Grassmann variables that still have to be integrated over. Writing the decoration factor in 
terms of the ip variables it takes the form 

8\{l^) - {fa)). (A.64) 

The decoration factors can then go on take part in further Grassmann integrations. 
This brings us to the final operation, 

C. MHV doubly glued to MHV: After several applications of operations A and B 
one will be left with a loop that connects two MHV vertices with some extra dependence 
on the integration variables in factors of the form 5 4 (...). 

Let us consider the example 

J d 4 r) h d 4 r) h 5 8 (q Kl + q K . 2 + X h r) h + A ?2 ^ 2 )5 4 ((/^) - (lip k ))5 8 (-X h r) h - X h r] h + q Ka ) 

= 5 8 (q) J d\ l( i\ 2 5\m - (l^))5\-X hVh - X hVl2 + q Ka ). (A.65) 

Here the ip variables are defined so that ipi —ip^ = qx 2 +^i 7 7h = — Qkx — \ 2 Vh- If we further 
define ip^ and ipi 2 via 

ip h - ip h = q Ks , ipi - iph = qx 2 (A.66) 

then we have that 

(im = (iM, m t ) = m 2 ). (a.67) 
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Hence we can reexpress if){ in terms if the 'external' if) variables, 



{hj) h )\i 2 - (jgjh)Aji 
(kh) 



(A.68) 



The final factor in (1A.65j) then just contributes (hh) 4 and we obtain 



5 8 {q)8 4 ((lif) l )(l 1 l 2 ) + M h )(l 3 l) + Mb) (Hi))- 



(A.69) 



Thus we see that the decoration factor develops more terms in the delta function. 

This calculation easily generalises to more complicated decoration factors. For example 
let us imagine that we have already performed one loop of Grassmann integrations. Then 
we will obtain a decoration factor of the type in (1A.69j) . When this appears under another 
integration of type C we obtain 



d A r] h d A r] h 5 8 (q Kl + q K2 + q Kz + \ h r) h + \ h r] h ) 

S 4 {(la^l a }(hk} + (lbA)(lila) + (lM(Ub))$ 8 (-KVh - KVh + QK t ). (A.70) 



5 8 (q)S 4 ({l a lP la )(l b h)(kk) + (hAWaHhh) + (lM(lih)(lah) + MMWJb)). (A.71) 



Generically there will be several decoration factors all depending on the internal if) 
variable. In this case all of them develop more terms under the Grassmann integration. 

Once all Grassmann integrations have been performed one can replace the external if) 
variables with dual 6 variables. The resulting expression is then clearly dual conformal 
covariant because every dual 9 is only ever contracted with the neighbouring A from the 
original diagram and A variables are only ever contracted with neighbouring A variables. 

The remaining step is simply to count the dual conformal weights of every factor we 
have produced throughout the calculation. Doing do one finds that they combine with 
the weights of the bosonic factors in precisely the right way. This also shows that all 
contributions have the same weight and so the result for the leading singularity is dual 
conformal covariant. 

B Two-loop integrals 

B.l Twistor identity to relate numerator factors 

In section [3] we gave an example for an identity that can be used to express pentabox and 
double pentagon integrals in terms of our preferred basis. The identity shown there is a 
special case a more general identity, 



(ll + l) = ajx(j-l j) + a j2 (J j + 1) + a kl (k - 1 k) +a k2 (k k + 1) +a jk (jk) + a^jk) , (B.72) 



The integral can be performed just as before using ( 1A.68j) as we obtain 



where 




(B.73) 
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The expansion coefficients are 
where 

a n = (jj + lkk + l)(k - 1 kll + 1) - (j j + 1 k - lk)(kk + lll + 1) , 

a j2 = -(j - ljkk + l){k - lkll + 1) + {j - Ijk - 1 k)(k k + 1 1 1 + 1) , 

a kl = (j-ljll + l)(jj + lkk + l)-(j-ljkk + l)(jj + 111 + 1), 

a k2 = -{j-ljll + l){jj + lk-lk) + {j -Ijk -lk){jj + 111 + 1), 

a jk = -{j - lj j + 11 + l){k - 1 k k + 1 1) + {j - 1 jj + 1 l){k - 1 k k + 1 1 + 1) , 

a^=-{jkll + l). (B.75) 

This can be used to simplify any of the pentaboxes or double pentagons appearing in any 
two-loop amplitude. 



"xy 



(j -ijj + ik)(jk-ikk+iy 



(B.74) 



B.2 Definition of the two- loop integrals 

Here we define the two-loop integrals appearing in (13.291) . The double pentagon integral 
with "magic numerator" is given by 



h;2;ml ~ I 2 s 2 



dZ^BdZcD 



(13 AB) (46O0) (6134) (1245) (2356) 
(61 AB) (12 AB) (23AB) (34AB) X 

1 



x- 



(ABCD) (34CD) (AhCD) (56CD) (61CD) 
and the remaining integrals are defined as 



(B.76) 



and 



f dZ AB dZ CD (i7i 2 )~ 2 (36CD)(6123)(1234)(1256) 

6;2;m 2 J (QiAB) (12 AB) {23 AB) (ABC D) {23C D) {34C D) {5QC D) (61C D) ' [ ' ' 

f dZ AB dZ CD {iix 2 Y 2 (35CD) (6123) (1234) (1245) 

6;2 ;m3 J (QiAB) {12 AB} (23 AB} (AB C D) (23C D) (34C D) (45C D) (5QC D) ' [ ' } 



f dZ AB dZ CD (ivr 2 )- 2 (1245) (6123) (3456) 

6;2;7 J (Q1AB} (12 AB} (23 AB}(ABCD)(3ACD}(45CD}(5QCD} ' 1 ' ' 

f dZ AB dZ CD {m 2 y 2 (1234)(6123)(2345) 

6;2;15 J (Q1AB} (12 AB} (23 AB} (ABC D)(23CD}(3ACD}(A5CD}' [ ' ' 
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B.3 Results for generic kinematics 

In the small m 2 expansion we obtain the following results. 



h:2;7 



6;2;15 



o(i), 

— — log 4 m 2 + log 3 m 2 lo£ 



2 2 2 

x 13 x 24 x 26 
x 14 x 36 



^ log 3 ™ 2 log'" ' 



„2 



rf 2 ,-y»2 

x 35 x 26 
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Here we used integrals like 



L 



dzx z Y I (-z)Y I (l + z) 



Re(*)=-l/2 



1 — X 



log(x) 



(B.86) 



We also computed the log 2 (m 2 ) analytically. One can check that the infrared divergent 
terms in ( 14. 35ft are correctly reproduced. 



B.4 Analytic results in the limit u — > 

Here we give the analytical results for all integrals for x 2 i+2 = 1 ,x 2 i+3 = u~ x l 2 as u — > 
(after having taken the regulator limit m 2 — > 0). We denote U = logw and L = logm 2 . We 
obtain 

hmi = \ uA + u2 Y + 2U( 3 + ^n 4 + 0(m 2 ,u), (B.87) 
hm2 = -^+ 1 -L 3 U+ 1 -L 2 (-U 2 -7r 2 ) + 1 -L(U 3 + 3U7r 2 -12C 3 ) 

- f/ 2 ^ - 2U( 3 - ^ + 0(m 2 , u) , (B.88) 

+ (^ 4 + ^U 2 + 4&U + ^) + C(m 2 , u) , (B.89) 
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and 



h = \l* + L*U- 2 ^L> + l(- 1 -UZ- 2 -£u + 4( 3 

+ (j- 6 u4 + 2 ~y u2 - 2U ^ + 70 + °^ 2 ^ > ( R9 °) 

^15 = y 2 L* - \l?U - ^L 2 + \L (f/ 3 + 4Ltt 2 - 48C 3 ) 

—f 4 - y ^ + 4f/Cs - |j) + 0(» 2 , «) • (B.91) 



Putting everything together, we obtain 



Mf = ^L 4 + L 2 ( 9 -U 2 3Cs£ + f ^I/ 4 + ^U 2 + ^ ) + 0[m\ u) , (B.92) 



and hence 

Jim Tlf («, «, «) = y log 2 « + ^ + 0(u) , (B.93) 

which is quoted in the main text. 

We have similar results for the limit u — > oo. 

C Box integrals for one-loop MHV amplitudes 

In order to write down the one-loop n-point MHV amplitude up to 0(m 2 ), we require the 
so-called "one-mass" and "two-mass easy" box integrals (in our case with uniform internal 
mass m). They are readily evaluated and are given by 

**V«, 4, x\ A ) = log 2 £ + log 2 £ - log 2 £ - log 2 X ^-%- 

x 46 x 15 x 14 J/ 46 ° 

9 \ / 2 \ 

X 14 \ ot; h x 14 \ , /n /2\ 



x 15 J \ X 46 



2Li 2 1 - ^ - 2Li 2 1 - ^ + £>(W) . (C.94) 



and 



f 2me (xL4,4 5 4) = -lo^^-io^^ + io^^ + io^^ + iog 2 ^ 

^14 ^36 ^13 ^46 -^36 
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^f) -2Li 2 (l-^ 
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+2Li 2 ( 1 - - 2Li 2 1 - ) + G(m 2 ) . (C.95) 



Here F lm = x| 6 a; 2 5 / lm and F 2me = {x 2 13 x 2 A6 - x\ A x\§)I 2me . Note that both functions can be 
identified with their dimensional regularisation counterparts (up to a constant in the finite 
part of F lm ) when expanding the latter in e and setting the dimensional regularisation scale 



ji = m 2 . 
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